1 


FAILURE  DISTRIBUTIONS  OF  SHOCK  MODELS 


GARY  GOTTLIEB 


TECHNICAL  REPORT  NO.  181 
October  18,  1976 


SUPPORTED  BY  THE  ARMY  AND  NAVY  

UNDER  CONTRACT  N00014-75-C-0561  (NR-042-002)  V’H' 
WITH  THE  OFFICE  OF  NAVAL  RESEARCH  ^ 


Gerald  J.  Lieberman,  Project  Director 


’<j  DEC  15 

LTL'lSlbli 

0/^  A 


Reproduction  in  Whole  or  in  Part  is  Permitted 
for  any  Purpose  of  the  United  States  Government 

Approved  for  public  release;  distribution  unlimited. 


DEPARTMENT  OF  OPERATIONS  RESEARCH 
AND 

DEPARTMENT  OF  STATISTICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA 


^ “ I • . ,1 


iwm— .i<pi  n m-  <■ 


0,  Non- Technical  Summary 

In  this  paper^  a single  device  shock  model  is  studied.  The 
model  we  consider  consists  of  a single  device  which  experiences 
shocks  which  come  from  the  outside  environment.  An  example  is  a 
sensitive  electrical  component  which  occasionally  experiences  a 
large  electrical  surge  due  to  a malfunction  in  the  electrical  system. 
Each  of  these  shocks  can  render  the  device  inoperable^  and  will  at 
least  make  it  more  likely  to  fail  when  the  next  shock  arrives. 

It  is  often  important  to  classify  the  distribution  of  the 
time  to  failure  for  some  item.  What  we  attempt  to  do  in  this  paper 
is  to  examine  the  one-device  shock  model  and  to  consider  the  condi- 
tions on  the  shock  process  and  on  the  ability  of  the  device  to  with- 
stand shocks  which  allow  us  to  do  such  classifying.  We  try  to  find 
conditions  on  the  process  which  allow  us  to  place  the  t ime- to- fai lure 
distribution  of  the  model  into  one  of  the  common  reliability  theory 
classifications. 
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1.  Introduction  and  Summar; 


In  this  paper^  we  shall  study  a simple  one-device  shock  model  and 
attempt  to  establish  various  qualitative  features  of  its  time- to- failure 
distribution. 

Let  L be  the  time  of  failure  of  the  device  and  let  H(  t)  = P(L  > t)^ 

and  H( t)  = 1 - H( t) . We  suppose  that  the  device  is  subject  to  a series 

of  shocks  whose  arrival  process  is  a renewal  process  with  interarrival 

distribution  F.  We  further  suppose  that  the  device  survives  the  first 

k shocks  with  probability  We  assume  that  ^ ^ ^ ^*2  — ’ ’ " ' 

Then  H(  t)  = y P [F^^^(t)  - F^  t) ) ] ^ where  F^^^  is  the  kth 

k=0  ^ 

convolution  of  F. 

The  general  approach  of  this  paper  is  to  make  various  assumptions 
about  the  sequence  and  the  sequence  of  functions  {F^  ( t) -F^  t) } 

or  fj  the  density  of  F,  and  to  then  try  to  classify  H( t) ^ the  time  to 
failure  distribution.  The  classifications  common  in  reliability  theory 
include  Increasing  Failure  Rate  ( IFR)  distributions,  Increasing  Failure 
Rate  on  the  Average  ( IFRA)  distributions  and  distributions  with  Polya- 
Frequency  of  order  2 ( PFg)  densities. 

The  general  setup  of  this  paper  largely  follows  that  of  [ 1]  and 
[3].  In  [3],  the  shock  renewal  process  was  taken  to  be  Poisson,  that 
is,  F was  assumed  to  be  exponential;  and  in  [l],  the  shock  process  was 
assumed  to  be  non-homogeneous  Poisson,  In  addition,  we  will  apply  some 
of  our  results  to  the  stochastic  wear  process  introduced  by  Morey  in 
[6],  A stochastic  wear  process  is  an  increasing,  non-negative  stochastic 
process,  where  F(t,x)  = < x)  is  Totally  Positive  of  order  2 (see 

Section  3).  We  will  now  discuss  some  of  the  results  of  [3]  and  of  our 
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results.  If  the  reader  is  unfamiliar  with  some  of  the  terms^  he  can 


refer  to  Sections  2 and  5. 

BeloWj  we  note  the  major  results  of  [jl,  which  can  be  extended^ 

in  some  natural  fashion^  to  the  case  of  a more  general  renewal  process. 

Many  of  the  other  results^  which  deal  with  distribution  classifications 

such  as  New  Better  than  Used  (NBU)  and  New  Better  than  Used  in  Expectation 

(NBUE)^  seem  to  be  valid  only  for  Poisson  renewal  processes. 

00  , 

Suppose  that  H(  t)  = 5L  ® (At)  /k.'  P where  1 = ^ — 

k=0  ^ 0 i ^ 

then 


i) 

If 

e.  = P, /P,  , is  decreasing 

in 

k = 1,  2,  ..., 

i.e.,  if 

k > 0}  is  a PF^  sequence. 

then  H is  IFR. 

ii) 

If 

r,l/k  . j , 

P|^  IS  decreasing  in  k = 

1, 

2,  . . . , then  H 

is  IFRA. 

iii) 

If 

Pk  = Vl  - \ ""  P^2» 

the  density  of 

II  is  PF^ 

We  will  extend  each  of  the  above  results  to  a more  general  renewal 
process . 


For  i)  ^ we  assume  that  F^'^^(t)  - F^"^^^^(t)  is  TP  for  ii)  ^ 


assume  that  f^'^^(t)  is  TP^  and  in  iii),  we  assume  that  f^'^^(t)  is 


TP, 


For  i)  and  ii),  no  further  assumptions  are  made  on  the  {Pj^} 
sequence^  but  for  iii),  we  assume  that  ~ 1 decreasing. 

In  i),  ii)  and  iii),  we  assume  that  7.  ff’^'^^(t)  - ( t)  )|'^ 


n=0 


is  respectively  log  concave  for  | £ [O, l],  IFRA  for  | £ [0,1]  and 

has  a log  concave  density  for  | £ [O, l]. 

In  each  case  we  get  the  same  result  as  in  [3]  and  we  also  show 

that  the  assumptions  on  the  quantity,  7 - F^ ( tj^*^ 

n=0 


are 
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necessary.  We  also  find  a large  class  of  densities  which  satisfy  the 


various  assumptions. 


For  the  stochastic  wear  process,  we  study  the  hitting  time  of 

to  some  random  barrier  X.  We  find  necessary  conditions  on  the 

process  Z and  on  X so  that  the  distribution  of  T„  where 

= inf[t  > o|z^  > X),  is  either  IFR,  IFRA  or  has  a PF2  density 

Zt 

Suppose  Z^  is  a stochastic  wear  process  with  E[|  ] log 

concave  in  t for  | € [O, 1].  Then  if  F( t,x)  = P(Z^  < x)  is  TP^ 

and  X has  a PFg  density,  or  if  f(t,x),  the  density  of  F(t,x),  is 

TP^  and  right  continuous  in  t and  X has  an  IFR  distribution,  we 

show  that  has  an  IFR  distribution. 

Also,  if  Z^  is  strictly  increasing  and  f( t,x)  is  TP^  and 

right  continuous  and  X has  a PF2,  decreasing  density,  we  show  that 

Tj,  has  a PF^  density.  In  all  of  the  above  cases,  we  show  the  necessity 

Zt 

of  the  condition  that  E[|  ] be  log  concave. 

Zh 

In  addition,  if  E[|  ] is  IFRA  for  | £ [O, l],  if  F(t,x)  is 

TP^  and  X has  a PF^  density,  we  show  that  is  IFRA. 


2,  Distribution  Classes  in  Reliability  Theory 

In  this  section,  we  define  some  of  the  classes  of  distributions 
used  in  reliability  theory.  Throughout,  F and  H will  be  distributions 
on  IR^  with  F(0)  =0.  f and  h will  be  respectively  the  density 
of  F and  H (when  the  distributions  are  absolutely  continuous).  F 
and  H are,  respectively,  1-F  and  1-H.  F^"^^  is  the  nth  fold  con- 
volution of  F and  f ^ the  nth  fold  convolution  of  f. 

will  be  a decreasing  sequence  of  non-negative  numbers  with  P^  = 1. 


When  appropriate^  we  will  have 


I 


H(t)  = I (F^")(t)  - P . 

n=0 

F is  said  to  be  IFR  if  F(  x+t)/jF(  t)  is  decreasing  (when  defined) 
in  t > 0 for  all  x > 0. 

If  F(  t)  ^ is  decreasing  in  t > 0,  F is  said  to  be  IFRA. 

It  is  well  known  that  F IFR  implies  that  F is  IFRA, 

F is  discrete  IFRA  if  the  support  of  F is  the  non-negative 
integers  and  F(n)^^'^  is  decreasing  in  n (n  an  integer), 

F is  discrete  IFR  if  the  support  of  F is  the  non-negative 
integers  and  (F(n))  is  PFg  (see  next  section)  on  the  integers,  A 
sequence  fa  ),  n > 0 is  discrete  PF_  if  and  only  if  a ,/a  is 

decreasing  in  n. 

If  F is  absolutely  continuous,  F IFR  is  equivalent  to 
r(  t)  = f(t)/F(t)  increasing  in  t > 0 ( F(  t)  > 0)  for  some  version  of 
f. 

5,  Summary  of  Total  Positivity 

Consider  any  function  K : X x Y ->]R,  where  X,  Yc  IR,  Define 

the  function  K^p^(x,y)  = detHK(x^,  p on  the  set 

j=i;.’!!> 

AP(X)  yAP(Y)  where  AP(X)  is  the  set  of  all  P-tuples 

for  which  X^  < i = 1,  P-1  and  X^  € X.  Define  AP(Y) 

similarly, 

*■ 

Throughout^  increasing  (decreasing)  means  non-decreasing  ( non- increasing) , 


I 

1 


I 
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Definition.  K(x,y)  is  (strictly)  totally  positive  of  order  r TP^ 
(STP^)  if  Kjpj(i,y)  > 0 (>  0),  P < r,  V(  x,y)  £ AP(X)  x AP(  Y)  , If 
K(x^y)  is  TP^  (STP^)  of  all  orders^  it  is  said  to  be  TP  (STP  ). 


Definition.  If  K(x^y)  is  TP^  and  admits  a representation  K{x^y) 
= P(y-x)j  we  say  that  0 is  a Polya  Frequency  Function  of  order  r 

( . 


A function  is  PF2  if  and  only  if  it  is  log  concave,  F is  IFR 
if  and  only  if  F is  PF2,  If  F has  a density  f which  is  PFg,  then 
F is  IFR. 


We  now  note  some  other  characterizations  of  PF^  functions  and 


IFRA  distributions.  From  [2]^  F is  IFRA  if  and  only  if  for  each 
A > Oj  F( t)  - e has  at  most  one  change  of  sign^  and  if  one  change 
of  sign  actually  occurs^  it  occurs  from  + to 

From  [3],  f 0 PF^  on  IR^  if  and  only  if  for  all  ^,6, 


- 9t 

f( t)  - ae  changes  sign  at  most  twice^  and  if  two  sign  changes  occur^ 


1 


negative  to  positive  to  negative 


Examples 


We  now  introduce  the  notion  of  variation  diminishing.  Let  F( t) 
be  defined  on  I where  I is  an  ordered  set  of  the  real  line.  For 
our  purposes,  I will  always  be  E''', 

Let 

S"(F)  = sup  S"[F(t^),  F(tJ]  , 

where  the  sup  is  taken  over  t,  <t_<"**<t,t.el  m arbitrary 
and  S Xj^]  is  the  number  of  sign  changes  of  the  sequence, 

zeros  being  ignored. 

Let  K(x,y)  defined  on  X x Y be  Borel  measurable.  Assume 

that  / K(x,y)  du(y)  exists  for  every  x in  X.  u is  taken  to  be 
Y 

a sigma- finite  regular  measure  on  Y.  Let  f be  bounded  and  Borel 
measurable  on  Y and  let 


g(x)  = / K(x,y)  f(y)  u(  dy)  , 

Y 

The  proofs  of  the  following  three  theorems  can  be  found  in  [4], 

Theorem  1.  If  K is  TP^  and  the  above  conditions  are  fulfilled, 
then  S ( g)  < S ( f)  provided  that  S'(f)  < r-1. 

Also,  if  f is  piecewise  continuous,  then  f and  g exhibit 
the  same  sequence  of  signs  when  their  respective  arguments  traverse 
the  domain  of  definition  from  left  to  right. 
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Theorem  2.  If  K(x,y)  is  TP^  on  (O,oo]  x (0,«>]  and  integrable  over 
y with  respect  to  a sigma- finite  measure  on  [O^oo]^  then  the 

iterates^  L(n^y)  are  TP^  for  n = 2^  3^  y ^ where 


L(l,y)  = K(0,y) 

00 

L(n,y)  = / L(n-1^  |)  K(|,y)  dn(|)  , n = 2,3,4,  ... 

0 


Theorem  3.  If  K(|,Y)  is  TP^  on  X x Y and  L(t,n)  is  TP^  on  Y x Z, 
where  X,  Y and  Z are  subsets  of  the  real  line  and  a is  a sigma 
finite  measure  on  K,  then 


M(^,n)  = / K(|,t)  L(T,n)  da(  t)  , | G X,  n £ Z 

Y 

is  TP  . / \ assuming  the  integrals  converge  absolutely. 

min(  m,  r j 

Lemma  1 . If  f is  a PF^  density  concentrated  on  [0,oo],  and  if  F 
is  the  associated  distribution  and  F^*^^  the  nth  fold  convolution  of 
F,  then  F^  ( t)  - F^"'^^^(t)  is  TP^.  f^"^(t)  is  also  TP^. 

Proof.  Apply  Theorem  2 with  K(x,y)  = f(y-x),  d^(^)  = d^  to  see  that 
f^")(t)  is  TP^. 

Apply  llieorem  3 with  K(|,t)  = f(T-|),  L(x,m)  = 
da( x)  = dx  to  see  that  1 - F( t)  is  PF2. 

Finally,  note  that  the  product  of  two  TP2  functions  is  TPg  and 
then  apply  Theorem  3 with  K(|,x)  = (1  - F(|-x))l^^  and  L(x,n)  = f^”^(x) 
and  da(  x)  = dx. 
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4 . The  Shock  Renewal  Process 

Let  us  more  formally  consider  the  shock  renewal  process.  In 
continuous  time^  let  be  the  number  of  shocks  which  have  arrived 

by  time  t^  and  if  we  assume  that  shocks  only  come  at  integer  times^ 
let  be  the  number  of  shocks  which  have  arrived  by  time  n. 

Then  for  ^ £ ( 0, 1) ^ 


E[|  "]  = E P(N  = n) 

n=0 


E - F^"+^)(t)) 


n=0 


= 1 + E F^"^(t)  , 


n=l 


E[^  ] = and 


lim  = 1 + lim  E,  F^")(t) 

^ t ->  00  n=  1 


t 


,00  / , (by  bounded 

= 1 + E F ”'^(t)  convergence) 

n=  1 t — > 00 


..  1 

F.  5"  . 

- n=l 


Nt,  Nf 

For  1=0,  let  E[|  ] = P(N^  = O) , Note  that  E[|  ] is  continuous 
in  Our  major  assumption  through  the  remainder  of  the  paper  is 

^2) 


Nt  Nn 

that  E[^  ] (E[|  J)  is  either  log  concave  ( PFp)  or  IFRA  (discrete 


IFRA). 
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We  list  some  examples  below: 


(1)  If  F(t)  = 1 - e"'^  , then 


= 7^  (.M)  ^ g-^t(l-5) 

n=0  "• 


(2)  Let  F(  t)  be  a gamma  distribution  with  rate  1 and  paraite  ter  2. 

2 

Let  f be  the  density  of  F.  Let  p = 

wAi  7 p(n)(,)  p2n  ^ 

P n=l 


Taking  Laplace  transforms  with 


g(s)  = / e'®*^  E[|  *^]  dt 


8(s)  = - + / e-®"  7 f('^)(t)  dt 

P 0 n=l 


I P^-1  1 ” -St  ...  2 ^ ^n 

P n=l  0 
P^ 

2 ? 

- i P i ( l+sl 
s 2 s ■ 2 

1 - — 


-n 


1 , (p2.i)  1 ^1 

(1+s)^  - p'^J 


Reinvertingj 
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E[|^']  = ^ [(UP)  ^ (^.1)  . 

A simple  calculation  will  show  that  this  function  is  log  concave, 

(5)  Consider  a renewal  process  on  the  non-negative  integers 
with  interval  distribution  where  F(  n)  = 1 - p*^^  0<p<l  q=l-p 

Such  a process  is  Markovian.  It  can  be  shown  that  N = J’  S where 

k=0  ^ 

are  independent^  identically  distributed  with  E(Sj^  = 0)  = q 
= 1)  = P.  Therefore^ 

k=0 

= q"( W P^/q)"  . 

Nn 

Clearly^  E[|  ] is  a PF^  function  on  the  non-negative  integers. 

(Ij.)  Let  G(  n)  = F'^Vn;^  where  F is  as  in  (5),  If  n is 

n 

the  discrete  renewal  process  associated  with  then 

[n/2j 

N = 7,  s 

k^O 

where  {S^}  is  as  in  (5).  Letting  p = we  have 
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E[| 


,n,  , 2k  (n-2k) 

>-^  (2J  (PP)  * 

k=0 


1 / n N / „^2k+l  n-(2k+l) 


(pp^-q)'"  ^ (pp-q)""  . ( 1 ^ i) 


This  sequence  is  PF^, 

Section  7 contains  a discussion  on  how  to  construct  distributions 
Nt 

for  which  E[|  ] is  log  concave. 

Theorem  4.  Let  F have  a density  f 
for  all  discrete  IFRA  sequences  {Ej^} 
decreasing  in  t for  all  | £ [O, l]. 

Also,  F^")(t)  - F^"'^^)(t)  TP^ 
f be  PF2  . 

Theorem  5 . Let  F be  a distribution  function  whose  support  is  the 
positive  integers  and  suppose  that  f(n)  = F( n)  - F(n-l)  is  PF^. 

Then  H(  n)  is  discrete  IFRA  for  all  discrete  IFRA  sequences 

Nn  ^ 

if  and  only  if  E[|  ] ' is  decreasing  in  n for  all  | £ [O, 1], 

Theorem  6.  Suppose  that  F^'^^(t)  - F^'^'''^^(t)  is  TP^.  Then  H(  t) 

Ni- 

ls IFR  for  all  PF^  sequences  only  if  E[|  ] is  log 

concave  ( PFg)  for  all  | £ [O, 1], 


which  is  PFg.  Then  H is  IFRA 
Nt  1/t 

if  and  only  if  E[|  ] ' is 

can  replace  the  condition  that 
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Theorem  7 . Let  F be  a distribution  function  whose  support  is  the 
positive  integers  and  suppose  that  F^*'^(n)  - is  TP^. 

Then  Hfn)  is  discrete  IFR  for  all  PF^  sequences  (P  } if  and  only 
if  E[|  "]  is  PFg  for  all  4 € [0,1]. 


Theorem  8.  Let  Pj^  = k > 1.  Suppose  that  f^"^(t)  is  TP^. 

TheOj  H(  t)  has  a PFg  density  for  all  PFp,  decreasing  sequences  {p  ] 

if  and  only  if  g(t^5)  = E f is  log  concave  for  all 

n=  1 

[0,1]. 

Before  we  prove  the  above  theorems,  we  must  first  prove  a couple 
of  technical  lemmas.  Lemma  2 is  needed  for  the  proofs  of  Theorems  6 
and  7 3nd  Lemma  3 for  Theorem  8. 

First,  we  introduce  the  function, 


gU,t^,  tg)  = E[|  for  I (£  (0,1],  tg  > t^  > 0 


Let  k = inf{n  > 0 : P(N  = n)  > 0].  We  defir 


g(0,  t t ) = lim  g(|,  t t ) 

= P(N  = k )/P(N  = k ) . 

2 ‘'1  '"l  '^1 


Lerrana_2.  If  F^  ^(t)  - F^  ^(t)  = P(N^  = n)  is  TPg,  then  for  any 
® for  which  H(  t^^)  >0,  there  exists  a ^ £ [0,1]  for 
which  g(|,  t^,  tg)  = H(tg)/ii(t  ). 


Proof.  g(  1,  tg)  = 1 and  g is  continuous  in  | £ [O^l].  So^ 

if  we  can  prove  the  existence  of  a | £ [O,  1]  for  which 


0 < g(^,  tg)  < H(t2)/ii(tp  < 1 = g(l,  t^,  tg)  , 


then  by  the  continuity  of  g^  the  letnina  is  proved. 

Now,  g(0  t t ) = P(N  = k )/P(N  = k ),  so  it  suffices 

^ Cg  Cl  Cl  Cl 

to  show  that  g(0,  ti,  t^)  < H( t2) /H( ti) , that  is, 


P(N  = k ) ii(t  ) < ii(t  ) P(N  = k ) , 
‘^2  1 1 1 


P(N^  = k^  ) H(ti)  = £ P(N^  = kj.  ) P(N^  = £)  P^ 


2 1 


£ > k "2  "1 

1 


< S P(N  = k ) P(N  = Z)  P 
£ > k *^11  2 * 


— t. 


P(N  = k ) H(t  ) . 
1 1 


The  inequality  follows  as  P(N^  = k)  is  TPg. 

The  following  lemma  is  proven  in  greater  generality  than  needed 
for  this  section.  The  more  general  form  is  needed  for  Section  8.  The 
reader  can  skip  this  proof  for  the  moment  and  just  refer  to  the  corollary. 

We  will  let  ^ be  some  o-finite  measure  on  IR^,  g some  non- 
negative, decreasing  function  on  IR^  and  f(t,x)  some  function  on 
IR^  xIR^  which  is  TPg  on  IR^  x A where  A is  the  support  of  n,  and 
which  is  right  continuous  in  x for  each  t. 
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Lemma_^.  For  f,  g,  ^ as  above,  if  / f( n( dx)  is  decreasing 

0 


in  t,  then  for  t^^  < t^. 


f(t2,x)  n(dx)  / fCtg^x)  g(x)  ji(dx) 


I i 0 


w OO 

I f(tpX)  ii(dx)  / f(t^,x)  g(x)  u(dx) 


< lim 

I t 1 


/ f(t_,x)  ^x(dx) 


/ |i(dx) 


where  we  assume  that  all  denominations  are  non- zero  and  all  integrals 


finite. 


Proof.  Let 


V = U(x,y)  = f(t  x)  f(t  ,y)  - f(  t y)  f(  t x) 


We  will  first  prove  the  right-hand  side  inequality.  As 


lim  / f(t,x)  f/  ^,(dx)  = / f(t,x)  n(dx)  , 
f f 1 0 0 


we  must  show  that 


CCtgjX)  g(x)  ji(dx)  / f(t.,x)  dx) 


</  M(dx)  / f(t  x)  g(x)  ^(dx) 

0 0 


'■Mi'll  11" 


Now. 


0 > II  U(x,y)  (g(x) 

X > y > 0 


- 8(y))  l^(dx)  V(dy) 


I I U(x,y)  g(x)  n(dx)  V(dy) 
X > y > 0 


+ II  U(x,y)  g(x)  p.(dx)  V(dy)  . 

y > X > 0 


Substituting  for  U(x^y)  and  separating  each  of  the  above  integrals 
gives  us  the  desired  inequality. 

Now^  we  consider  the  left-hand  side  of  the  inequality.  Let 


k(t)  = infix  > oil  f(t,x)  ^(ds)  >0]  . 

0 


If  M.(  (k(  tj^)  ])  > 0,  the  inequality  follows  inmediately  from  the  fact 
that  f(t,x)  is  TP^  and  g( x)  >0.  Suppose  that  n(  {k(  tj^)  ))  = 0, 


lim 

I ‘ 0 


r » 

/ f(t_^x)  |i(dx) 

0 


I f( t,,x)  n(dx) 
LO 


= iim 
I I 0 


I n(dx)] 

k(ti) 


,I  i 


c-k(  ti) 


< lim  lim 

e i 0 5 i 0 


^k(  t^)+e 

I ^(dx) 

k(t^) 

k(  t ) +e  7^  V 

/ f(t  x)  ^(dx) 

1-  / ^ V A 


>(  t^) 

f(  tg,  k( t p + e) 


< lim  ,,  , , — r 5- 
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p(dx) 


r 


-«  "II  M\i\^nvmmmit9mmmm 


As  f is  TP, 


2^ 


i I 


; 


f(  k(  tj^)+e)  k(tj^)+e 


g(x)  ^i(dx) 


< 


I 

k(tp 


+e 


g(x)  ^(dx) 


Letting  6 -»0  on  both  sides,  and  noting  the  above  string  of  inequalities 
shows  that 


00 

^ / f(to,x)  g(x)  ji(dx) 

1 1 m — — ■ — ^ 0 

— • 00 

f i n(dx)  / f(t  x)  g(x)  n(dx) 

U Q i 


lim 
i i 0 


as  required. 

We  will  now  prove  the  needed  corollary.  Let  be  a decreasing 

non-negative  sequence.  Let  f( t)  be  the  density  of  some  non-negative 
random  variable. 


i 

I 


I*  £<">(.)  is  iP^,  d.„s.si„s, 

there  exists  a ^ £ [o, 1]  such  that 


h(  t. 


h(  t ) _ 
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Proof.  We  can  apply  Lemrr.a  3 with  g( k)  = p^^,  f(t,n)  = and  u 

counting  measure  on  the  positive  integers.  We  can  extend  g to 
and  f to  R xlR  in  an  appropriate  manner.  By  Lemma  3^ 


k=l 


<:^<  u„ 


on  ' V.l'  (-  V '' 

k=  1 k=  1 ^ 


00 

The  corollary  follows  by  the  continuity  of  f ^ ( t)  in  |. 

Ifcl 

Proof  of  Theorem  U. 


H(t)  = ^ [F^")(t)  - P 

n=0 


H(t)  - £[5^*^]  = E [F^")(t)  - F^"+^)(t)]  (P  - 

n=0  n ^ • 


,Vn 


As  P^  decreases  by  assumption,  the  sign  of  p - changes 


2' 


sign  at  most  once,  and  if  so,  from  + to  -.  As  f is  PF 
F ^(t)  - F^  ^ ( t)  is  TPg  by  Lemma  1.  So,  by  Theorem  1,  H(  t)  - 
has  at  most  one  sign  change  and  if  it  does  change  sign,  it  is  from 
+ to  - . 

We  suppose  that  H( t)  = 1 - H(t)  is  not  IFRA.  It  then  follows 
from  [2  ],  p.  89^  that  there  exists  a 0 < < 1 and  a t^  for 

which  H( t)  crosses  from  below  at  t = t^.  Now 


N. 

lim  E[|  O]  ^ = 0)  < H(  t ) 

, i 0 '^0  - 0 
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As  E[^  is  continuous  in  5 and  E[1  ] = 1 >H(t)^  it  follows 

N t 

that  there  exists  a £ [O,  l]  for  which  ^ H(  t^)  = 


- V Nt 

Therefore^  H(  t)  crosses  1 from  below  at  t = t . So 

Nt  ° 

the  sign  of  H(t)  - E[|^  ] goes  from  - to  +.  This  is  a contradiction 

so  H(  t)  must  be  IFRA. 

We  are  given  then  H(  t)  is  IFRA  for  all  sequences  {Pj^)  where 
1 = Pq  - ^1  - decreasing.  Let  = 1*^,  0 < ^ < 1, 

H(t)  = E (F^")(t)  - F^'^+^^Ct))  . 

n=0 

Nt  Nt  i/t 

So  1 - E[^  ] is  IFRA  (E[|  ] ' decreasing)  as  H is  IFRA. 


Proof  of  Theorem  . 

^ f(n)  = F(n+1)  - F(n)  is  PF^  by  assumption.  By  Lemma 

f^‘^)(n)  is  TP^, 


H(n)  = E P(N^  = *^)  I’k  = \ » 


k=0 


k=0 


N. 


H(n)  - E[|  “]  = E (P  - . 

k=0  ^ 


?k  ■ ^ changes  sign  at  most  once  and  in  the  order  + to  There- 

in 

fore,  by  Theorem  1,  H(n)  - E[^  ] changes  sign  at  most  once  and  in 

the  order  + to 

Suppose  H(n)  is  not  discrete  IFRA.  Then  we  can  find  a 
^0  ^ *^0  ^ ° which  H(n^)  < and  H(n^+1)  > 

We  can,  as  in  the  previous  proof,  assert  that  there  exists  a | € [O  l] 

Hn  N_  ^ ' 

for  which  IqO  > e[^  0]  > H( n^) . 

As  and  is  decreasing,  E[^"0+1]  < 

r Nm  . ° ' 

J goes  below  H(n)  at  n = n +1.  That  is,  the  sign  of 

H(n)  - E[|  ] changes  from  + to  - as  n goes  from  n to  n +1 

0 0 ’ 

^ Let  P^  = I and  proceed  as  in  Theorem  5. 


Proof  of  Theorem  6. 


00  . 

H(t)  = E [F^"\t)  - F^"+^)(t)]  P 
n=0 

H(t)  - aE[|”'^J=  E [F^")(t)  - F("+^)(t)]  (P  - a?")  . 

n=0 

As  {P^j  is  PFg,  is  decreasing,  so  the  sign  of  P^  - a|" 

changes  at  most  twice  and  in  the  order  - f As  F^'^^(t)  - F^”'''^\t) 
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Nt 

is  TP^,  we  can  assert  by  Theorem  1 that  the  sign  of  H( t)  - aE[|  ] 
changes  sign  at  most  twice  and  in  the  order  - + 

We  know  that  H is  IFR  if  and  only  if  for  any  a and  the 
number  of  sign  changes  of  H(  t)  - number  at  most  two^  and  if  there 

are  two,  they  are  in  the  order  _ + 

Let  us  suppose  the  contrary.  That  is,  assume  that  there  exists 
an  a,  | for  which  the  sign  of  H(  t)  - a|*"  goes  from  - to  + to 
Refer  to  the  below  diagram.  By  Lemma  2,  we  can  find  a b,  5 with 

r ^*^1.  . Nt  to  n: 

bE[|^  n = ai  ^ = H(t^),  bE[|^  '^  = H(  t2) . As  bE[|  "]  is 

Nt  t 

log  concave  by  assumption,  we  have  that  for  t^  < t < tg,  bE[|^"]  > a|\ 

It  is  clear  that  there  exists  t . t,  . t . t < t,  < t,  < t_  < t for 

a»  b'  c^  a 1 b 2 c 

which  H(tJ  - bE[|^'^^]  >0,  H(  t^)  <bE[|^‘^‘’],  and  H(  t^)  >bE[|j^^‘=], 

This,  however,  is  in  contradiction  to  the  initial  result  about  the 

Nt 

sign  change  pattern  of  H( t)  - aE[|  ] for  all  a,|.  So  H is  IFR. 

Again  let  P = |*^  and  the  converse  is  immediate. 

' n 
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t I 


Proof  of  Theorem  7. 

<=  A review  of  Lemma  2 shows  that  it  is  still  applicable  for 

discrete  time^ 

H(n)  = Z = Z P^^  , 


k=0 


k=0 


H(n)  - aE[|  = Z (P  - a?'")  . 

k=0 


As  Pj^  - a|  changes  sign  at  most  twice^  and  in  the  order  - + 

f k'l  . N 

and  as  f'  '’(n)  is  TP^,  then  by  Theorem  1,  H( n)  - aE[|  changes 

sign  at  most  twice  and  in  the  order  - i-  - . 

We  suppose  that  H( n)  is  not  PF^^  or  equivalently  that  there 

exists  constants  a',  1^,  n^,  n^,  where  n^  < n^  < n^  are  integers 

for  which  H(n^)  > a H(n2)  < a'  H(n^)  > a' 

Making  use  of  the  fact  that  any  two  exponential  curves  cross 

at  most  once^  we  can  find  a new  pair  of  constants^  a and  | for  which 


n2 


H(ni)  = a^^  , n(n^)  < a|^  , H(n,)  = a^^^. 

As  in  the  previous  theorem,  using  Lemma  2,  we  find  a b,  ^ 

N _ N 

with  bE[^j^  ^1  = H(n^),  bE[  5 j = H(  n^)  . We  proceed  as  in  the  pre- 
vious theorem  to  get  a contradiction.  So  H(n)  is  PF2. 

Let  P = 1". 
n 


1 


“>  / V 

h(t)  = z , 

n=l 

h(t)  - ag(t,|)  = Z (P^  - 

n=  1 


We  now  have  the  same  setup  as  in  Theorem  6.  We  now  may  use  the 
Corollary  to  Lemma  3 in  the  identical  way  that  we  used  Lemma  2 in 
the  proof  of  Theorem  6. 

^ Let 


Nt 

7.  Densities  for  Which  E[|  ] is  log  Concave 

As  was  mentioned  earlier^  the  question  of  finding  distributions 

with  the  properties  required  in  Theorems  1+  through  8 is  formidable. 

Nt 

Below,  we  identify  a class  of  densities  for  which  E[|  ] is  log 

concave. 


Definition.  Let  be  the  set  of  functions  Q on  [0,w]  with  the 

5s  k °° 

representation  Q(  s)  = ( e ® s^  II  (1  + A.s),  A.  > 0,  Z A.  < oo, 

n=l  ^ ^ i=l  ^ 

5 > 0,  k a non-negative  integer  and  Q(0)  = C, 


Theorem  9.  h( u)  is  a PF^  function  on  [0,oo]  if  the  reciprocal  of 
the  Laplace  transform  of  h is  in  ty 


Proof.  The  proof  of  this  theorem  can  be  found  in  [4]. 


Note  that 


= II  (A.+s)/nA, 
1 1 ^ 


n ( 1 + i/A.  S) 

1 ^ 


is  the  reciprocal  of  the  Laplace  transform  of  the  convolution  of  n 
exponential  densities  with  rates  A^. 

Let  G(t,^)  = 1 - 


Theorem  10.  Let  F have  a density  f*G(tj|)  has  a PF  density 

q(tj5)^  for  |£  [O, ij  if  the  Laplace  transform  of  f has  a reciprocal 
n n n n 

of  the  form  H (A.+s)/n  A.  and  ( II  (A.+s)  - II  A.)  has  all  real 

1 ^ 1 i=l  ^ i=l  ^ 

roots. 


Proof. 


q(t,e)  = E 5" 

^ n=l 


q^(s,|)  = I e'®  q(t,0  dt 

0 


( t\  _ IzL 

" I i-Us)i 


We  want  to  show  that  q( t |)  is  a PF  density,  and  as  it  is  clearly 
a density^  it  suffices  by  Theorem  9 to  show  that  q^^(s,|)  is  in 
for  all  I 6 [0, 1), 


(q^^s)  - I) 


= (n(Vs)  - ^ II  A^)  ( E A^)-^  . 

i=  1 i=  1 


Now,  (s,g)  is  in  if  p^(s)  = n(A^+s)  - | n A^  has  all 

real,  non-positive  roots.  For  | £ [0,l],  0^(s)  >0  for  s > 0. 

So  /i^(s)  has  no  positive  roots.  An  inspection  of  the  polynomial 
will  show  that  f)^  has  all  real  roots  for  | £ [O,  l]  if  and  only  if 
has  all  real  roots. 


Corollary . Let  f be  the  convolution  of  two  exponential  densities. 
Then  q(t,|)  is  PF^  for  | £ [0,1). 

=T^  (A2+S)  - 5 A^  A2))(A^A2)‘^ 

We  need  to  see  if  0^(  s)  = ( A^+s)  ( A2+S)  - A^^A2  has  all  real 
roots  0j^(s)  =■  s(s  r A^  + , so  the  conditions  of  Theorem  10  are 

satisfied . 


No^:  The  result  in  Theorem  10  only  holds  for  | £ [0,1),  but  it  can 

easily  be  extended  to  the  case  of  ^ = 1 by  continuity  or  by  noting 
that  q( t,5)  =0  for  | = 1. 


8.  Stochastic  Wear  Processes 


In  [6],  Morey  introduces  a class  of  non-decreasing  stochastic 
processes,  which  he  called  stochastic  wear  processes  ( SWP) . t > 0) 

is  a SWP  if  P(Z^  < x)  is  TP^  in  t and  x,  is  increasing  and 

Z^  is  non-negative. 

Morey  shows  that  if  Z^  is  a SWP  with  stationary  independent 
increments,  then  the  random  variable  T defined  as  infft  > olz  > xj 
is  an  IFR  random  variable.  He  further  asserts  that  for  any  non-negative 
random  variable  X,  is  also  IFR.  The  latter  result  is  not  true. 

Simply  take  to  be  identically  t.  note  that  P(Z  < x)  = I,  ^ , 

which  is  TPg  and  therefore  Z^  is  a SWP  with  stationary,  independent 
increments.  Also,  P(  < t)  = P(Z^  >X)  = P(X  < t) , so  IFR  implies 

that  X is  IFR,  which  is  clearly  not  true  in  general. 

From  now  on,  we  will  be  considering  a non-decreasing  stochastic 
process  Z^  with  = 0,  F( t,x)  = P(Z^  < x)  and  f( t,x)  the  density 
of  F for  each  t (if  it  exists). 

We  will  prove  the  analogues  of  Theorems  4 through  8 for  processes 
which  can  take  values  on  all  of 

First,  we  need  the  following  lemma  which  is  similar  to  Lemmas  2 

and  3. 


Lemma  4.  Suppose  that  F( t,x) 
Then,  for  all  t^^  < t^. 


is  TP^,  q > 0 and  p(  t) 


/ 

x= 


F(t,x)  g(x)dx. 


lim 
% i 0 


/ F(  t„,x)  dx 


x=0 


/ t,,x)  dx 

Ux=0 


Pi  tg) 

< VT  — v'  < lim 
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/ F( tp,x)  dx 


x=0 


/ F(  t x)  l’'  dx 
x=0 


Proof.  The  left-hand  side  of  the  inequality  follows  by  Lemma  ). 


Now,  F(t2,x)  <F(t^,x).  Also 


OC  00 


/ F(t,x)  I dx  = / / f(t,u)  du  dx 

x=0  x=0  u=0 


00  00 


/ / f(t,u)  du  dx 

u=0  x=u 


So. 


/ / f(u,  t)  1^^  dy  du 

u=0  y=0 


= -1/ln  i E[|  "]  . 


lim 

I t 1 


/ F(t2,x)  dx 


x=0 


/ F(t,,x)  dx 


l-x=0 


. u.  oiAtii-L  , 1 

i t 1 -l/ln  5 E[| 


Let  G(x)  = P(X  <x),  G(x)  = 1 - G(x).  Let  g(x)  be  the 

density  of  G.  Let  H(  t)  = P( T > t) , H( t)  = 1 - H( t) . Let  h( t) 

be  the  density  of  H( t)  (if  it  exists).  Then, 

00 

i)  H(t)  = P(T  > t)  = / F(t,x)  g(x)  dx. 

x=0 

ii)  H(t)  = P(T  > t)  = f f(t,x)  G(x)  dx. 

x=0 

If  is  strictly  increasing  in  t, 

00 

iii)  h(  t)  = / f(  t,x)  g(x)  dx. 
x=0 


2? 


Theorem  11.  If  F( C^x)  is  TP^,  then  H( t ; is  PF^  (H  is  IK  fur 
all  g PFg  if  and  only  if  E[|  1 = - / f(t,x)  f dx  'In  t is 

X 0 

log  concave  in  t for  all  ^ £ [0,1  . 


Proof. 


H(t)  - a / F(c,x)  f,’'  dx  / Ffi,x)  ;g(x)  - dx  . 

x=0  x^O 


H(t)  - a F:'t,x)  dx  changes  sign  at  most  twice  and  in  tlie 
x=0 

order  - + - as  g( x)  is  PF^  and  F(t,x)  is  TP  ^ Theorem  1 . Now, 


apply  Lemma  * and  argue  as  in  Theorem  5. 


Let  g(x)  = I , 


00 

Note  that  E[|  ] is  log  concave  if  and  only  if  T F(  t,x) 


is  log  concave. 


We  state  the  next  three  theorems  without  proof.  They  each 
depend  either  on  Lemma  j5  or  Lemma  4. 


Theorem  12.  Let  f(t,x;  be  right  continuous  for  each  t and  TP^.  Then, 
H(t;  is  log  concave  (H  is  IFR)  for  all  G(  t)  PF.^  if  and  only  if 

Zt 

E[|  ] is  log  concave  for  all  | £ [0,l], 


Theorem  Ij.  Suppose  that  is  strictly  increasing.  Let  f(t,x) 

be  right  continuous  for  each  t and  TP^.  Then,  h( t)  is  PF^  for 
all  g(  t)  PF^,  and  decreasing,  if  and  only  if  E[|  ] is  log 

concave  for  all  | £ [O, Ij. 
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Theofeni  1^.  Let  f( t,x)  be  right  continuous  and  TP^.  Then  H(  t) 
is  I FRA  for  all  G(  t)  IFRA  if  and  only  if  E[^  Vt  decreasing. 

Example.  Let  be  a non-homogenous  Poisson  Process  with  E(N^)  = A(  t) 

Consider  the  compound  renewal  process  whose  arrival  process  is  and 

whose  jump  size  distribution  is  M with  M'O)  . 0.  Assume  that  M 
has  d density  m.  Let  be  the  size  of  the  nth  jump.  Let 


if  N|.  > 1 


if  Nj.  = 0 


P(Z^<x)=  Z M^'^^fx) 


Let  F(t^x)  _ P(Z^  < x)  . Let  f(t^x)  be  the  density  of  F on  (O^oo], 
We  have 


, 


X > 0 


/ \ k 

Suppose  m(x)  is  PF  . Then  m^ ( x)  is  TP,.  Since  *^) 

> 3 k.' 

is  TP^,  we  know  by  Theorem  3 that  f( t,x)  is  TP^  on  t > 0,  x > 0. 
Extend  f(  t,x)  to  t > 0,  x > 0 so  that  it  is  still  TP  Now, 


E[f/'^J  = „here 


m = / m(x)  dx  . 
^ x=0 


So  E[f,  I is  log  concave  if  A(  t)  is  convex. 


Therefore,  by  Theorem  12,  if  A is  convex  and  m PF^,  the 
hitting  time  distribution  of  to  any  random  barrier  with  an  IFR 

distribution  is  IFR. 


9.  Markov  Processes 

The  questions  we  have  been  asking  about  hitting  time  distri- 
butions have  very  simple  answers  if  the  process  is  Markovian  and 
spatially  homogeneous.  All  of  the  below  is  taken  directly  from  [5J. 

Henceforth,  is  a Markov  Process,  which  is  spatially  homogeneous. 

We  let  p^(s,x)  be  the  density  of  P(Z^  < x|z^  = 0). 


Theorem  1^ . If  0^(s,x)  is  TP^  in  s and  x for  each  t and  if 

00 

f(x)  is  PF  then  c( t)  = / t x)  ffx)  dx  is  PF„. 

r\  6. 


Proof. 


c(  t+s)  t ) s,  x)  f(x)  dx 


00  X 
0 0 


/ I f(u+0  du  d^ 

0 0 


As  f is  ^('-'+0  is  RRg  (see  [kj).  By  a slight  extension  of 

00 

Theorem  ^ (see  [ ]) , / 0^(s,u)  f(u+^)  du  is  RRg,  so  as  ^^(t,^) 
is  TPg,  by  Theorem  5,  c( t+s)  is  RR^.  But  c( t+s)  RR^  implies  that 


c( t)  is  Using  Theorem  I5,  we  can  get  the  same  results  as  in 

Section  8 for  temporally  homogeneous  Markov  Processes  with  fewer 


and  weaker  assumptions.  Note  that  for  appropriate  choices  of  f,  c 
can  be  interpreted  as  the  hitting  time  distribution  or  density  of  the 
Markov  process  to  a random  barrier. 
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In  this  paper j'^w^considerj  a single  device  shock  model.  The  device 
experiences  shocks  from  the  environment, ^^®ch  of  which  can  render  the  device 
Inoperable. "we  find  sets  of"  conditions,  oinj  the  shock  process  and  S^the 
ability  of  the  device  to  survive  shocks/so  that  the  time  to  failure  distri- 
bution of  the  device  falls  into  one  of  the  common  reliability  theory  classi- 
fications. He  extend  these  results„to  the  case  where  the  shock  process  can 
be  viewed  as  continuous.  •''c  a-*  . .f<  .< 
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